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ABSTRACT 

Cluster  formation  is  simulated  numerically  with  discrete  velocity  Boltzmann  model 
in  two  space  dimensions.  The  model  exhibits  cluster  coagulation,  fragmentation,  and 
transport.  It  evolves  on  two  different  scales  obtained  from  an  elastic  and  inelastic 
collision  Knudsen  numbers  s  and  n  respectively.  For  flow  impinging  on  a  wall  with 
specularly  reflective  boundary  condition  these  scales  appear  both  analytically  and  nu¬ 
merically. 
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§1  Introduction 


In  recent  paper  Slemrod,  Grinfeld,  Qi  and  Stewart  [8]  have  given  a  model  and  one  dimensional 
numerical  simulation  for  a  discrete  velocity  gas  exhibiting  coagulation  and  fragmentation. 
For  example  such  a  model  might  illustrate  same  of  the  grosser  features  of  nucleating  droplets 
that  make  up  fog  or  an  aerosol. 

The  model  is  based  on  the  rules  of  discrete  velocity  kinetic  theory  and  many  of  the  ideas 
using  in  developing  the  model  are  to  be  found  in  an  earlier  model  of  Monaco  and  Pandolfi 
Bianchi  [2,  3].  The  new  feature  of  the  model  of  [8]  is  that  it  allows  clusters  of  particles  with 
unlimited  cluster  size  while  the  model  [2]  did  not. 

The  model  of  [S]  contains  two  length  scales  based  on  an  elastic  Knudsen  number  £  and 
inelastic  Knudsen  number  /t.  Since  inelastic  collision  processes  proceed  much  slower  than 
elastic  ones  we  expect  0  <  e  <<  fi.  This  in  turn  implies  at  least  for  steady  flow  that  we 
should  see  different  regimes  of  gas  flow  in  a  neighborhood  of  a  boundary. 

In  this  paper  we  consider  this  issue  analytically  and  numerically.  Indeed  in  a  sequence  of 
photos  of  computer  generated  data  we  have  able  to  see  the  boundary  layers  rather  sharply 
completely  consistent  with  the  analytical  analysis. 

§2  The  discrete  velocity  coagulation-fragmentation  model 

We  consider  a  discrete  velocity  gas  of  identical  particles  each  of  mass  m  contained  in  volume 
V  C  R^.  A  point  in  V  is  identified  with  its  Euclidean  coordinates  {x,y,z).  The  particles 
are  grouped  into  clusters  possessing  1,2,  ...,M  particles  where  the  cluster  of  size  M  is  the 
largest  cluster  that  can  be  crammed  into  V.  The  clusters  move  in  V  with  fixed  momenta 
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=  m(c,0,0),  £2  =  m(-c,0,0),  £3  =  m(0,c,0),  £,  =  m(0,-c,0),  £5  =  m(0,0,c), 

£5  =  m(0,0,-c),  £0  =  m(0,0,0).  A  cluster  made  up  of  a-particles  will  be  called  an 

P 

a-cluster.  It  is  clear  that  an  a-cluster  with  momentum  £j  has  velocity  u"  = 

As  the  clusters  move  in  V  they  collide  in  a  binary  fashion.  Both  elastic  and  inelastic 

collisions  are  allowed.  In  terms  of  the  momenta  £  the  elastic  collisions  will  be  represented  by 
list  £i  +£2  =  £3  +£4  =  £5  +  £5  the  inelastic  collisions  are  defined  by  £  +£2  =  £0, 

£3  +  £4  =  £0.  £  +  £  =  £0,  £.  +  £0  =  £.  (*  =  •••’  6)-  Notice  elastic  collisions  conserve 

mass,  momentum,  and  energy;  inelastic  collisions  conserve  mass  and  momentum  but  not 
energy.  We  denote  by  n^(x,y^z,t)  the  number  density  of  a-clusters  with  momentum  £  at 
a  point  {x,y,z)  €  V  at  time  t  >  0,  i.e.  the  number  of  clusters  in  this  class  per  unit  volume. 

A  collision  of  an  a-cluster  with  momentum  £  and  a  ^-cluster  with  momentum  £•  which 
yields  a  ^-cluster  with  momentum  £.  and  a  7-cluster  with  momentum  £  will  be  represented 
by  (n",nj)  — >  (nf,n^).  This  notation  allows  us  to  write  the  allowable  elastic  collisions  as 
follows. 

1.  Mechanical  collisions:  (”35^4) 

with  similar  statements  for  (713,714 ),  {n°,n%). 

2.  Exchange  collisions:  (a)  “head-on”  collisions:  {n°,n2)  (^4)^2)  (”^>^2) 

(pro6.|),  (773,774)  {prob.\),  (773, 77J)  {•pTob.\),  {n^,nQ)  {prob.\),  (775,77^)  {prob.\) 

with  similar  statements  for  (773,774)  and  (771,775).  (b)  “angle”  collisions:  (77i',77f)  —i’ 
(77“,  773)(pro6.|),  (774,773)  {prob.^)  with  similar  statements  for  (77“, nf)  and  (773,775). 

We  only  allow  inelastic  collisions  of  the  Becker-Doring  type  [4,  5,  7],  i.e.  where  an 
a-cluster  may  gain  or  lose  a  1-cluster  in  coagulation  or  fragmentation  respectively.  The 


coagulation  of  a  1-cluster  with  momentum  P_j  to  form  an  a  -F  1-cluster  with  momentum 
Pjt  is  represented  as  (nl,n“)  while  the  fragmentation  of  an  a  -|-  1-cluster  with 

momentum  Pj^  into  a  1-cluster  with  momentum  Pj  and  an  a-cluster  with  momentum  Pj 
will  be  denoted  by  (n^+^)  (nl,np.  With  this  notation  the  allowable  Becker-Doring 
inelastic  collisions  are  as  follows. 


1.  “Head  on”  coagulation:  (n},nf  ^),  (nf  \n^),  (n^,<  ^),  (n| 
(n^^nJ) (n“). 


2.  “Moving  cluster  coagulates  with  rest  cluster”:  (nj,no“'),  (77."  \77q)  — >  77“,  j  = 

1  6 


3.  Reversalofl  (fragmentation):  «) ->  (n}, 7iJ“^)(pro6.|),  tz^)  (pro6.|), 

(pro6.J),  (na'^TzJ)  (pro6.|),  (pro6.|),  {nt~\nl)  {prob.\)  if  a  >  2; 

(no) (n},7zj)  (pro6.|),  (n^,7z*),  (pro6.|),  (72^,71^)  (pro6.i). 

4.  Reversal  of  2  (fragmentation):  (tz^)  (7zj,7Zo“^)  {prob.\),  (7zJ,7z““*)  {prob.\),  j  = 

1,2,..., 6  if  Q  >  2;  {n])  -7  (nj,7zi),;  =  1,2,..., 6. 

The  rate  equations  governing  the  motion  of  clusters  are  given  by  the  transport  equations 


dn°  c  dn° 

dt  ^  a  dx 
dn^  c  dn2 

dt  a  dx 

dn%  c  dn^ 

dn^  c dn° 

dt  a  dy 

dn%  c  dn° 


= 

=  + 

=  + 

=  E°  +  I^, 
=  Et  +  I!, 


l<a<M 


(2.1) 
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driQ  c  dn% 
dt  a  dz 

dt 

The  calculation  of  EJ  has  been  given  in  [2,  3]  according  to  the  rules  of  discrete  velocity 
kinetic  theory,  i.e.  these  terms  are  proportional  to  collisional  cross  sectional  areas;  the 
relative  velocity  of  the  particles  before  collision;  the  probability  of  each  admissible  collision; 
both  number  densities  of  the  colliding  particles.  Clearly  Eq  =  0, a  =  We  record 

the  rest  of  Ej  and  the  inelastic  collision  terms  /“  in  Appendix. 

It  is  a  straightforward  exercise  to  show 


E^  +  Et  +  E^  +  E^  +  E"  +  =  0  ,  a  =  1, ...,  M  , 


M 


M 


£r  -  £?  =  E  ^3°  -  =  H  Et  -  =  0 , 

a=l  Or  =  l  0=1 

M 

E  a(/o”  +  It  +  4“  +  4“  + 1°  +  4“  +  4”)  =  0 . 

0=1 

M  M  M 

E  /f  -  4“  =  E  4”  -  4“  =  E  4“  -  4“  =  0  ■ 


a=\ 


0=1 


a=l 


Next  define  the  quantity: 


Ja  =  ”o(“o.*i'*l  + 

+(1  -  +  a^',1n2  +  +  «o>5  +  «o,i<) 

+  (1  —  d"  ‘^0,1^3  ^4  d"  ^o’.l’^S^e) 


d-  <jto’,l”4’^3  d- 


tl  „o+l  t2  „o+l  43  „o+l  _  lA 

^0+1^1  ^0+1^2  ^a+l^3  ^0+1^4 


0+1 


®o+l“5  ®o+l”6 


for  1  <  a  <  M  —  1  , 


Jm  =  0 . 


4 


It  then  follows  that 


IS  +  /f  +  IS  +  It  +  It  +  It  =  ^-1  1 


Let  us  define 


n“  =  n“7l 


where  7^  is  a  typical  value  of  the  density  Yl%oTla=i  ^  ~  ^I'l  V  ~  vl'i  ^  ~  ^  ~ 

where  L  is  a  typical  macroscopic  length.  Hence  n",  5,  y,  z  are  dimensionless  quantities. 
Then  for  example  in  equation  (2.1)  we  see 

dn^  Idn^  _  L 

dt  a  dx  cR  ^  ^ 

Hence  we  see 

dh°  I  dn°  _  r  2^  ( pa)  ,  I'  ra 

where  we  define  as  with  c  =  1,  n®  replaced  by  n“,  and  ri  =  1  in  the  definition  of  cr^/s. 
Next  set  e  =  Ij Lr\'R.  which  is  the  dimensionless  elastic  Knudsen  number. 

/“  possesses  both  coagulation  and  fragmentation  coefficients.  It  is  natural  to  assume 
the  coagulation  coefficients  scale  in  a  similar  way  to  the  elastic  collision  coefficients,  i.e. 
they  are  proportional  to  cr\.  By  consistency  we  then  assume  the  fragmentation  coefficients 
are  proportional  to  ct\%. 

Thus  setting 


2-hk  ^ 

= 

r 

=  cr?TC6i+,4 
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where  fi  is  the  dimensionless  inelastic  Knudsen  number,  and  is  If  with  nj  replaced  by 
nj,  a{’*,  6^+1  replaced  by  the  above  definitions,  we  see 

^  +  iM  =  lL  +  £. 

dt  a  dx  £  n 

Finally  drop  the  over  bars  and  we  have  derived  system 


1  dnf 
a  dx 
1  dnf 
a  dx 
1  dnf 
a  dy 
1  dn'^ 
a  dy 

a  dz 
1  dnf 


^  +  (a) 

£  fl 

rpQ  ja 

^  +  (6) 


t  n 

pot  ja 

e  fi 

K  +  l 

£  ft 

-I-  1± 

£  H 


I  <  a  <  M 


We  note  that  since  in  our  model  the  total  number  of  particles  is  always  conserved,  a  cube 
of  volume  will  have  =  Mof  where  A/q  is  the  number  of  particles  at  f  =  0.  Of  course 
AIq  must  be  less  than  or  equal  to  M.  Hence  we  have  £  =  L^/v^Mq.  In  the  Boltzmann  limit 
[1]  for  rarefied  gases  we  take  Mq  —*  oo,  >  0,  vol(F)  — <•  oo  with  Mo  finite.  This 

motivates  our  study  of  the  case  of  M  =  oo,  £  finite.  However  in  the  Boltzmann  limit  for 
dense  gases  we  take  Mq  —>■  oo,  j  0,  with  Mo  finite  i.e.  since 

^  {njLfMo  ’ 


we  consider  £  — »  0,  the  fluid  dynamical  limit. 
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Next  record  in  addition  the  equations  for  transport  of  o-clusters 


+  |-(iV"u;“)  =  ,  2<a<M 


+  fziNW)  = 


where 


j=0 


N°w‘^ 


-«  -  <)  , 
a 


-(n?  - »;)  , 

a 

-«  -  O ; 


the  equation  for  conservation  of  mass 


where 


d  d  .  .  ^  r  \  ^  f  \  n 


=  2aiV“u“  =  X:K-n?), 

a=l  cr=l 

M  M 

=  £aiV“t;“  =  £K-n"), 


=  ZaNV  =  ZK-<)’ 
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and  the  equation  for  conservation  of  linear  momentum 


+  +  O  =  0. 

fi  f)  ^  1 

|(P-)  +  ||>r+<)  =  0. 

If  we  define  the  macroscopic  symmetric  tensor  11  by 

a 

Uyy  =  -/9u''  +  ^^(77"+<), 

n,-,  =  -/9ty^  +  -(n?  + '^e)  > 

Q 

Hi-y  =  —pUV,  fix;  =  -piiW,  Ilyj  =  -pVW 
then  the  conservation  of  linear  momentum  may  be  expressed  in  the  familiar  form 


(2.6) 


dt 

d 


d 


d 


d 


d 


dx 

d 


dy 

d 


d: 

d 


0  . 


In  this  paper,  we  consider  only  the  two  space  dimensions  by  setting  n'j{x,y^z,t)  = 
nj{x,y,t)^  making  the  symmetry  assumption  nf{x,y,t)  =  nQ{x,y,t)  and  forcing  u;  =  0.  As 
noted  above  in  the  Boltzmann  limit  we  take  M  =  oo. 


§3  A  simple  scaling  analysis 

Consider  the  transport  equation  (2.2)  in  the  Boltzmann  limit  M  —  oo.  We  wish  to  study 
the  behavior  in  the  neighborhood  of  wall  x  =  0.  First  we  rescale  x,  t  and  set  x  =  i  =  f 
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where  T  is  a  typical  dimensionless  time.  Substitution  of  this  scaling  into  (2.2a)  yields  for 


example 

1  dnX  1  dn^  ^  /f 

T  di  ect  dx  e  fi 

Hence  for  small  time  when  T  =  £  we  find  to  leading  order  the  system  is  governed  by  the 
elastic  collision  system 


dn^ 

di 

1  dn° 
a  dx 

=  E7, 

(a) 

dn^ 

di 

1  dn2 
a  dx 

=  ^2  » 

(b) 

dn^ 

di 

171  or 

—  ^3  ? 

(c) 

dn° 

di 

= 

id) 

Out 

di 

=  Et, 

(e) 

dn^ 

=  0, 

if) 

dt 


oo 


(3.1) 


as  long  as  y  is  outside  a  layer  of  width  e  from  the  y  =  0  wall. 

On  the  other  hand  for  larger  times  when  |:  <<  1  we  recover  the  equations 


1  dn^ 
a  dx 

=  Et, 

(«) 

1  dn% 
a  dx 

=  E^, 

ib) 

(3.2) 

fi  dn^ 

T  di 

= 

(c) 

0  =  £S 

=  Et 

=  Et. 

(3.3) 

For  intermediate  time  T  =  /i,  (3.2c)  yields  the  dynamic  equation 


f)n° 

{T  =  n) 


(3.4) 
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while  for  large  time  (3.2c)  T  »  ft  yields  the  constraint 


0  =  IS  {T»tt). 


(3.5) 


In  the  neighborhood  of  the  vertex  x  =  0,  y  =  0  we  rescale  y  in  a  similar  manner  y 


= 

ff’ 


T  =  €  yielding  the  system 


+ 


di 

dn^ 

di 

dn^i 

di 
di 


1 

Q  dx 
1  dnl 
a  dx 

a  dy 

iM 

a  dy 


=  Et,  (a) 

=  Et,  (6) 

=  £3  ,  (C)  1  <  Q  <  00 

=  E- .  (d) 

=  E°  ,  (e) 

=  0,  (/) 


(3.6) 


Again  for  larger  times  ^  <<  1  we  recover  the  equations 


I  dn"^ 
a  dx 
1  dnl 
a  dx 
1  dn° 
a  dy 

iM 

a  dy 


«r.  (“) 

£?,  ((.) 

,  (c)  1  <  a  <  00 

Et,  (d) 


0  =  £?.  (e) 


fi  driQ 

f~M' 


■'0  • 


(/) 


(3.7) 


Again  for  T  =  ^  and  T  »  j.i  we  obtain  (3.4),  (3.5)  respectively  from  (3.7f). 

Notice  of  course  for  y  independent  motion  (3.7)  reduces  to  (3.2),  (3.3).  Thus  for  large 
time  T  »  fi,  the  motion  in  an  £— neighborhood  of  the  wall  x  =  0  for  one  dimensional  flow 
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and  in  a  square  with  sides  of  length  e  near  i  =  0,  y  =  0  for  two  dimensional  flow  is  governed 
by  steady  elastic  collisions  combined  with  the  constraint  Iq  =  0. 

Now  we  must  do  the  analysis  on  the  larger  7^  scale.  The  crucial  step  is  to  rewrite  the 
cluster  balance  laws  (2.3),  (2.4)  in  the  scaled  variables  i  =  f,  x’  =  y*  =  In  this 
notation  (2.3),  (2.4)  become 


2  <  a  <  oo  ,  (3.8) 

(3.9) 

On  the  same  time  scale  T  =  s  the  elastic  terms  dominate  and 

we  find  from  (2.2)  the 

spatially  homogeneous  system 

dn° 

^  £"  =  0  ,  1  <  i  <  6,  1  <  a  <  oo 

dt  ^  ^ 

(3.10) 

=  0  ,  1  <  Of  <  oo  . 

dt 

(3.11) 

On  the  intermediate  scale  T  =  fx  spatially  inhomogeneities  enter  via  (3.8),  (3.9)  and  we 
obtain  the  equations 


|:(A^“)  +  ^(iV“«“)  +  ^(7V%“)  =  Ja-i-Ja,  2  <  a  <  oo  (3.12) 
dt  ox"  ay" 

+  +  (3.13) 

dt  dx-^  dy"  ^=1 

combined  with 

=  IS  (314) 

dt 

from  (2.2g).  In  addition  we  have  the  constraint  that  the  solution  evolve  on  an  elastic 
Maxwellian 

EJ  =  0  ,  1  <  ;■  <  6,  1  <  a  <  oo  ,  (3.15) 
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from  (2.2a-e). 


We  recall  from  [8]  that  (3.15)  implies 


(n^7^“,n“,<,0  =  (e‘=^e-‘=^e‘=^e-<=^l)C'^  •  (3.16) 

For  longer  times  T  >>  /x  (3.12),  (3.13),  (3.14)  reduce  to  their  time  independent  steady 
state  versions  obtained  by  setting  time  derivatives  equal  to  zero. 

For  motions  where  x  >>  /x,  y  >>  t  >>  1^  the  fluid  dynamic  limit  becomes  the 

dominant  system: 


where  n"  evolve  along  the 


dp  d  ,  ,  5  ,  , 

Tt  +  a?*'’”’ 


d  ,  ,  5  ^  1/  a  ^ 


=  0. 

=  0, 

=  0, 


manifold  of  elastic-inelastic  Maxwellians  (again  see  [8]) 


(3.17) 


given  by 


£“  =  /?  =  0  ,  0  <  i  <  6  ,  I  <  a  <  00  , 


Qa 

Qi 


(l,e=^e-"',e^^e 

rt0,-2 

b.---b2  ’ 

1, 


'^l)(Q„(ni)"  , 

a  >2 


(3.18) 


Cl?  C2,  tIq  are  functions  of  {x,y,t).  In  this  case  (3.17)  yields  a  system  of  three  conservation 
laws  in  the  dependent  variables  ci{x,y,t),  C2{x,y,t),  nl{x,y,t).  The  system  is  meaningful 
in  fluid  dynamic  limit  £  — >  0  for  finite  values  of  the  density  p,  i.e.  when  0  <  nj  <  and  2, 
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is  the  radius  of  convergence  for  the  power  series 


a=l 

Here  p,  =  'ET=i  is  the  saturation  density. 

Thus  we  conclude: 

1.  In  an  e-layer  near  the  wall  the  elastic  Boltzmann  like  regime  dominates  given  the 

(3.1)-(3.7). 

2.  In  a  wider  layer  of  width  /z  the  cluster  transport  is  taking  place.  This  layer  interpolates 

between  the  macroscopic  fluid  equations  (3.17)  and  the  microscopic  elastic  collision 
Boltzmann  like  dynamics  (3.1)-(3.7).  This  wider  ^  layer  is  governed  by  equations 
(3.8)-(3.16). 

There  is  one  additional  scaling  that  is  appealing  to  make  in  the  transition  layer  of  width 
-  pL  when  (3.12)-(3.16)  are  valid.  Specifically  we  follow  the  discussion  of  Zel’dovich  and 
Raizer  [10,  VI,  Section  18]  and  treat  a  as  a  continuous  variable,  1  <  a  <  oo.  This  yields  the 
rule  for  evolution  along  elastic  Maxwellians  as 

(ni(x,  j/,  a,  t),  n2{x,y,Q,  t),  y,oc,t),  n4(x,  y,  a,  t),  n5)(x,  y,  a,  t) 

=  (e'=‘,e-‘=Se‘=2,e“^Sl)C(x,y,a,t) ,  1  <  a  <  oo  ,  (3.19) 

where  Ci,  are  functions  of  (x,y,t).  Then  we  find 

6 

N{x,y,a,t)  =  Ylnj{x,y,a,t) 

3=0 

=  no{x,y,a,t) 

+C(x,  y,  a,  t){2  cosh  ci  +2  cosh  C2  +  2)  ,  (3.20) 
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(3.21) 


2 

N{x,y,a,t)u{x,y,a,t)  =  -(sinh(ci(x,j/,0)C'(a:,j/,Q:,<) , 

2 

N{x^y,Q,t)v{x,y,a,t)  =  —{sinh{c2{x,y,t))C{x,y,a,t),  (3.22) 

J{x,y,a,t)  =  a{a)no{x,y,a,t)C{x,y,a,t){2coshci+2cos\iC2  +  2) 
+a{a)nQ{x,y,  I,  t)C {x,  y,  a,  t){2  cosh Ci  +2coshc2  +  2) 
+3a{a)C{x,y,  l,t)C{x,y,a,t) 

+3a{a)C{x,y,  l,t)C{x,y,aJ)  (3.23) 

- b{Q)C {.i\y,  a  J){2  cosh  Cl  +2coshc2  +  2) 

—  (2 cosh  Cl  +  2cosh  c-2  +  2)  — (6(a)C(.r, y, a. <))  , 

Io{x,y,a,t)  =  3l){a)no{x,y,oi,t) +  3a{a)C{x^y,i,t)C{x,y,a,t) 

-3C{x,y,  l,t)—{a{a)C{x,y,a,t)) 

+3a{a)C{x,y,  l,t)C{x,y,a,t) 

-3C{x,y,  l,t)—{a{Q)C{x,y,a,t)) 

-a{a)C{x,  y,  1,  t)C{x,  y,  a,  0(2  cosh  Ci  +  2  cosh  cj  +  2) 
-|6(a)C(a:,y,a,0(2coshci  +  2coshc2  +  2)  (3.24) 

1  d 

-i(2coshci  +  2coshc2  +  2)—{b{a)C{x,y,a,t))  ,  a  >>  1  , 

Io{x,y,Ut)  =  -  r  ai0)no{x,y,l,t)C{x,yJ,t){2coshci+2coshc2  +  2) 

J0=i 

-i  {b{p)C{x,  y,  /?,  0(2  cosh  c,  +  2  cosh  cj  +  2) 

+-^  W)C{x,  y,  13, 0(2  cosh  Cl  +  2  cosh  C2  +  2)]  \  . 

op  ) 

we  have  conveniently  though  not  rigorously  regarded  1  as  a  small  parameter  when 


a  >>  1  and  used  the  approximation 


/(Qr  +  l)«/(a)  +  /'(a)  . 


If  we  combine  (3.23)  (3.24)  with  the  continuous  version  of  (3.12)  we  find  for  large  a  that  the 
governing  equations  become  the  conservation  laws 


d 


d 


dJ 


di 


no{x,y,a,t) 


Io{x,y,a,t) 


(3.25) 

(3.26) 


where  N,  Nu,  Nv,  J,  Iq  are  given  by  (3.20)-(3.24).  This  provides  a  system  of  equations  for 
the  variables  nl{x~,y', q.  i),  C’{x’,if, o,  f),  cj(.r'‘,  y’,  f),  C2(.t*,  ?/’,  i)  obtained  from  no,  C,  Ci, 
C2  by  the  indicated  change  of  variables.  When  combined  with  the  equations  of  conservation  of 
momentum  (2.6a, b)  written  along  (3.19)  we  have  four  equations  in  four  unknowns  depending 
on  the  four  independent  variables  i’,  ?/*,  a,  i.  Unfortunately  (3.25),  (3.26)  will  possess 
the  nonlocal  “boundary”  terms  C{x,y,l,t)  no{x,y,l,t)  which  must  be  obtained  from  the 
continuous  versions  of  (3.13),  (3.14): 


—.N{x,y,l,t)  +  —{N{x,y,l,t)u{x,y,l,t))  +  —{N{x,y,l,t)v{x,y,l,t)) 


/OO 

J{x,y,a,t)da  , 


—.no{x,y,l,t)  =  Io{x,y,l,t). 


(3.27) 

(3.28) 


A  continuous  approximation  of  (3.12)  would  envolve  solving  (2.5),  (2.6)  (subject  to  (3.19)), 
(3.25),  (3.26)  for  a  >  1  with  boundary  conditions  (3.27),  (3.28)  and  appropriate  initial 
conditions  for  n^,  C*,  cj,  c^. 
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§4  Numerical  simulations 


In  this  section,  we  present  the  numerical  simulations  of  a  gas  impinging  upon  wall(s)  in  a 
two  dimensional  square  box.  We  assume  specularly  reflecting  boundary  conditions.  Using 
the  discrete  velocity  model,  we  expect  to  observe  the  nucleation,  that  is  the  formation  of 
larger  clusters,  on  the  boundary  where  the  gas  impinges. 

The  numerical  implementation  of  the  model  is  straight  forward  except  for  the  method 
used  to  solve  the  partial  differential  equations  (2.2).  For  that  we  used  the  first  order  upwind 
scheme  to  discretize  (2.2).  Also  we  require  the  closure  hypothesis  Jm  =  0  to  enforce  mass 
conservation,  M  finite.  The  calculations  are  carried  out  at  U-W  Madison  on  an  HP710 
provided  by  NSF  under  a  SCREMS  grant.  The  color  pictures  were  created  on  an  SGI 
machine. 

The  kinetic  coefficients  in  the  inelastic  collision  process  are  computed  with  the  formulas 
used  by  Penrose,  Lebowitz,  Marro,  Kalos,  and  Tobochnik  [6]: 


For  the  case  a  =  1 


and 


2a^  =  -(874  +  1888a)^/^  ,  2  <  a  <  oo  , 

6 

(4.1) 

bo+i  =  2waaa  ,  2  <  a  , 

(4.2) 

(4.3) 

10,  =  0.010526  . 

(4.4) 

=  1(874  +  1888)^/^ 

6 

(4.5) 

^  1  o 

Wa=  Q  OC-l,Z. 

(4.6) 
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63  =  2w2a2  ,  h2  =  loifli  . 


(4.7) 


Qi  =  l  and  we  choose  Q2  =  exp(2.599),  Qz  =  exp(5.70S)  from  the  equilibrium  simulation  of 
Stauffer,  Coniglio,  and  Heerman  [9]  given  in  Table  1  of  [6]. 

The  initial  data  are  given  by  the  formulas 


(7z",n“,n“,n^,n“,0 

{no,n\,nl,nl,n\,nl) 


(0,0, 0,0,0)  ,  2  <  a  <  M 

(l,e^^e“'■^e‘=^e■'=^l)(5l7J‘  , 


(4.8) 

(4.9) 


Therefore  initially,  we  have  only  the  1-clusters  in  the  system. 

Choosing  Uq  =  1.0,  C\  =  2.0,  C2  =  0.0,  Knudsen  numbers  t  =  15  and  /.i  =  30  and  M  =  S, 
we  have  a  case  that  can  be  reduced  to  one  dimensional  model.  The  Fig.  l(a)(b)(c)(d)  are 
illustrated  the  density  distribution  of  the  gas  at  time  steps  50,  150,  250  and  350  respectively. 
In  all  of  gray  scale  photo  pictures,  the  variations  of  gray  scale  from  white  to  grey  to  black 
correspond  to  the  density  variations  from  zero  density  to  lower  density  to  highest  density. 
The  gas  impinged  upon  one  side  wall  and  was  drawn  off  from  the  opposite  wall.  The  Fig. 
2(a)(b)(c)(d)  shows  the  density  distribution  of  total  nucleated  clusters,  that  is  the  all  clusters 
except  of  the  1-clusters,  for  the  corresponding  steps  as  in  the  Fig.  l(a)(b)(c)(d).  The 
nucleation  process  are  clearly  visible  and  proportional  to  the  density  distribution  of  the  gas 
at  the  side  of  impinged  wall. 

The  Fig.  3(a)(b)(c){d)  show  the  density  distribution  at  steps  50,  150,  250  and  350 
respectively  for  the  gas  sending  to  a  corner  of  the  square  box.  Here  nj  =  1.0,  Ci  =  2.0, 
C2  =  2.0.  £  =  20,  /X  =  40  and  M  =  8.  The  Fig.  4(a)(b)(c)(d)  show  the  corresponding 
nucleation  density  distributions.  The  Fig  3  and  4  are  illustrated  the  similar  phenomenon  as 
in  the  Fig  1  and  2  except  that  in  the  Fig  3  and  4,  the  phenomenon  is  not  reducible  to  one 
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dimensional  model.  Both  computational  visualizations  show  the  formation  of  two  layers,  one 
nested  within  the  other  at  the  boundary. 

§5  Conclusion 

We  have  performed  a  numerical  simulation  for  the  discrete  velocity  kinetic  model  with  co¬ 
agulation  and  fragmentation  originally  described  in  [8].  The  computations  illustrate  the 
formation  of  a  boundary  layer  containing  the  large  clusters  (denser  vapor)  which  is  con¬ 
tained  within  a  thicker  layer  in  which  the  rarefied  vapor  is  being  converted  to  the  denser 
vapor.  Of  course  the  computations  only  reflect  the  two  length  scales  built  into  the  original 
model:  £,  //.  In  an  £— neighborhood  of  the  wall  the  elastic  collisions  dominate  and  coagula¬ 
tion  and  fragmentation  ceases  while  on  a  larger  neighborhood  and  away  from  the  wall  it  is 
the  inelastic  collisions  with  the  flow  evolving  along  the  manifold  of  elastic  Maxwellians  which 
dominate.  Since  it  is  the  inelastic  collisions  that  cause  the  nucleation  we  observe  nucleation 
occurring  through  the  wider  /i— layer.  Outside  the  wider  layer  the  macroscopic  conserva¬ 
tion  laws  of  mass  and  momentum  (isothermal,  compressible  Euler  equations)  are  valid.  This 
is  meaningful  in  the  infinite  limit  M  oo  as  long  as  0  <  p  <  p*,  pa  th®  saturation  density. 
In  the  p,  e  layers  we  have  no  such  restriction  and  p  can  exceed  p,. 
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APPENDIX 


Formulas  for  elastic  terms  EJ: 


—  (^3^4  +  HgTlg  2njn2} 

o  a 

6  0^ 

■i-n'^riQ  +  ti^tIq  —  5Uj  /^2 ) 


i)(nfn^  +  <”4  +«3< 


+  ^571^  +  TI^TIq  577j  n2  ) 

M 

+  Yi  C(X^0Rop{>An%  +  nfn"  +  nfnf  +  nfn^ 
/J/or 


"nf  -  <nj  -  <71?  -  n"nj)  , 


4  ^TofO^C  ^  ^ 
o  a 


+  E^(-  +  l)(- 

6  0’ 


i)(n"n^  +  n>J  +  njn“ 


+nfnf  +nfn|  -f>n%nl) 


M 

■  ^  C(T0,pRa(}{p-2 

0z^a 


n"  +  nfn"  +  nfnf  +  nfrig 


-njng  —  n^nf 


.^nj  -  njJi?  -  nlnl)  , 


4  ^00^ !  a„a  I 

- - (71,712  + 

o  a 

M _ 1 


«  -  27i3“n:: 


+  E  ^  (-  +  4)  ("3<  +  "?"?+ "f "? 

0^a  6  0' 

+77>?  + 71^77“  -57l>f) 

M 

+  Y1  c<Tc0Ra0{nin^  +  nln^  +  Tif  Tzf  +  tzJti" 

Pl^Ct 

-TlgTl,  -  773772  -  ^3^5  “  ^S^e)  > 

l^(nrn;  +  nX-2nX) 

o  a 

+  E  ^(-  +  4)(">J + <”? + "f"? 

/3^a  6  ^G  0^ 


ta 

+775775  +  77^776  -  577^774) 
M 

■  YL  C<^c.0Ra0{rAnX  +  77^77^- 

0^c 


nf  +  774775  +  7747 
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-n°n{  -  -  n“nf )  , 


Et 


^(TaaC 

3  a 

M 


(n“n2  +  ~  2^5 rig) 


+  E  ^(-  + 

6  'a 


+n>f +  nf<-5nX) 


A/ 


Pi^a 

-<nf  -  <nf  -  -  <nf )  , 


^6 


4  cf 


3  a 

M 


«nj  +  n?nj  -  2nfnS) 


+  L^(-  +  l)  (”>i  +  "?«?  +  "f "? 

rE  6  '^a  0' 


Pi^a 


+72371^  +  72372^  —  ^n^Tig) 


M 


+  X]  +  726  722  +  +  72.6  724 

P^ot 

-726  Tif  “  ^6^2  ■“  ^6^4  “  ^6^4) 


Here  CTap^  cRap  denote  the  collisional  cross  sectional  area  and  magnitude  of  relative  ve¬ 
locities  in  a  (nj,72^)  collision.  An  estimate  of  (Xap  is  easily  obtained  by  noting  that  it  is 
proportional  to  (r^  +  '^on'^p  radii  of  a-  and  ^-clusters  respectively.  Since  the  vol¬ 
ume  of  an  a-cluster  in  ^7rar?  the  radius  of  a  spherical  a-  cluster  is  and  hence  (Tap 

is  proportional  to  To  compute  Rap  we  note  that  the  velocity  of  an  rij 

cluster  is  ^  and  an  nf  cluster  is  Hence  Rap  =  (q  +  |)  “head  on”  collisions  and 
Rap  =  (^^  +  for  collisions  at  a  right  angle. 


Inelastic  terms  If  for  2  <  a  <  M: 


n  =  +  (1  -  '5»2)a!,’2i,inr'nJ  - 


+(1  -  -  ai>rn;  +  +  -fo-nK*') . 


^2 


—  ^  +  (^  “  ^02)00-1,1^2  ^”0  “  ^0^2 


L2„a 
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Oijj  0\  Q  COq 


+(1  -  S^M){-al'\n^n\  -  > 

=  <^a’-l,l’^3^0  ^  +  (1  ~  ^02)^0’- 1,1  ^3  ^”0  ~  ^a”3 

+(i  -  «»M)(-aj'»;  -  <-”"X + . 

=  ^a’-l,1^4’^0  ^  +  (1  —  ^o2)Oa-l,1^4  ^^0  ”  ^a^4 

+(1  -  «.«)(-4f,n:n5  -  <.i>;n;  +  56J+.T.J+'  +  i62„nS+')  , 
=  aa-l.l«5«r^  +  (1  -  ^a2)aa^,inr^nj  - 

+(i  -  -  4>?"5 + +  j<'°+i'‘r‘) . 

=  ®a-l,l”6'^o”^  +  (1  ~  ^a2)aa°i,ing  —  b^TlQ 

+  (1  —  ^QA/)(“<^aJ^6”5  “  ^a’,1^6^0  2^a+1^6^  "t"  g^Q+1^0^  )  ’ 

=  -6«n"  +  al\,n\nr"  +  4’^.in3<''  +  «"-i.i«5<"' 

+  (1  ~  ^a2)(Ol[o-l^?  ^^2  ®l!a-1^3  *^4  ^l,ot-1^5  ^e) 

+(1  —  <5oa/)(— a°’\nQn}  —  o-a^i^on^  — 

-«a,t»^0«4  -  «S”o^5  -  «a,l<^6  +  ^^a+l"l 

J-^A2  -0+1  I  ^;,3  „o+l  I  1^4  „a+l  ,  1l5  a+1 

+  9®«+1^2  +  2^a+1^3  +  ^  2^“+^”® 

4.1^6  „a+l\ 

+  2^Q+1^6  )  • 


Inelastic  terms  If  for  a  =  1: 


M—1  1 

t\  \~^  f  1.2  10.  1.0  1  0  /i  I  c  \/^l1  ,^1+)^  I 

A  =  -  2-/  \Ki3^1^2  +  -  (1  +  0/3i)(2^+/?^1  +  “^  °  n  ’ 

A/~l  1 

A'  =  -  E  -  (1 + . 

/?=1  ^  ^ 

M— 1  1 


M-1  1 

t\  r  4.3  1  j3  ,  4.0  1  /3  n  I  C  M  I 

A  =  -  (oi^n^n^  4-  afpn^n^  -  (1  +  hi){2^i+0^4  +  “^”0  ) /  > 
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M-1  1  /.O 

/  5,6  1  /?  I  5,0  1  /J  n  X  \f^k^  »,1+^_L 

=  ~  2^  |«i;;a”X  +  <^1,0^5^O  -  (1  +  ^PlK^^l+p^S  +  ';  ’ 

/?=1  ^  ^ 

1  Ifi  ^  ^ 

r  6,5  1  0  ,  6,0  1/?  /I  ic  ^/■^L6  t 

=  “  2^  (ai;^n6<  +  a{pne,n^  -  (1  +  hiR-^^i+p^e  +  "^"o  ))  , 

/?=1  ^  ^ 

M-l 

(  QA  I  0  .  0,2  1  0  .  0,3  1  0 

-  ~  1^  X^hP^O^l  +  ai!/3”0«2  +  Ol,/?«0% 

P=1 

+aj;jnjnj  +  ajjnjn?  +  alfpnln^ 

+^i+/3^4^^  +  ^i+;0”5^^  +  ^i+ZJ^e"*"^)}  • 


The  and  are  positive  kinetic  rate  coefficients  for  coagulation  and  fragmentation 
respectively.  The  satisfy  the  principle  of  detailed  balance  In  general 

^or’$3)  ^Q+i  should  depend  on  the  volume  V  but  in  the  classical  Becker-Doring  approach  we 
use  here  they  are  assumed  independent  of  V . 
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X  c^)  t-  tv 


-t-tro 


t-zro 


3  (i) 


-t 


±-  ITq 


'  c  ^  ^ 
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I.C) 


f  =  zro 


